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ABSTRACT 


Calculation  of  entropies  for  molecules  having  internal  rotations  are  in 
error  if  these  rotations  are  treated  as  vibrations.  Using  MNDO-calculated 
geometries,  entropies  of  a  set  of  molecules  were  calculated  assuming  all 
internal  motions  were  vibrations.  Errors  were  determined  for  these  molecules, 
with  particular  attention  paid  to  the  ones  containing  internal  rotations. 

Based  on  the  total  internal  partition  functions  calculated  in  this  manner 
compared  to  total  internal  partition  functions  calculated  by  assuming  a  free 
rotor  for  the  internal  rotation  we  developed  a  criterion  for  determining  when 
the  free  rotor  correction  needs  to  be  made  to  the  entropy  values.  The 
correction  greatly  improves  the  calculated  entropy  values  for  molecules  with 


internal  rotation 
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CHAPTER  1 
INTRODUCTION 


The  MNDO  (Modified  Neglect  of  Diatomic  Overlap)  computer  program  developed 
by  Profeasor  Michael  J.  S.  Dewar,  University  of  Texas  at  Auatin,  can  be  used 
to  calculate  vibrational  frequencies,  thermodynamic  quantities,  when  coupled 
with  appropriate  statistical  mechanics,  and  the  optimized  geometries  of 
molecules.  While  MNDO  calculations  are  fairly  representative  of  most 
experimental  entropy  (S°)  values,*  whenever  there  are  free  or  slightly 
hindered  internal  rotations  (such  as  a  methyl  rotor  adjacent  to  a  triple 
bond),  the  MNDO  program  produces^ slightly  erroneous  entropy  values  at  298.15  K 
fofvvmany  of  these  molecules.  Since  the  calculated  entropy  values  are 
important  for  determining  a  variety  of  other  values,  the  desire  for  a  higher 
degree  of  accuracy  is  obvious. 

A  considerable  amount  of  research  has  already  been  done  on  internal 

rotational  barriers  and  methods  of  modeling  them  in  various  molecules. 

Professor  Dewar  and  associates  have  written  many  articles  about  the  results 

obtained  from  MNDO  and  its  predecessor,  MINDO/3.  Some  of  these  results 
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include  the  computer-generated  optimized  geometry  of  molecules  and  entropy 
values  of  molecules  without  internal  rotations.*  Both  —  the**.  optimized 
geometry  of  molecules  and  the  entropy  values  for  the  molecules  without 
internal  rotations  are  very  close  to  experimental  values.*  However,  in 
molecules  with  free  or  only  partially  free  internal  rotations,  the  basic 
calculations  erroneously  model  these  internal  rotations  as  low  frequency 
vibrations.  These  low  frequency  vibrations  make  the  vibrational  partition 
function  contribution  much  larger  than  it  should  be,  which  subsequently 
results  in  an  incorrect  vibrational  entropy  contribution. 


This  paper  describes  one  method  of  correcting  MNDO  calculations  involving 
molecu?es  with  free  internal  rotations  (IR)  or  partially  hindered  rotations. 
It  consists  basically  of  removing  the  spurious  vibrational  mode's  (modes') 
contributions  to  the  vibrational  entropy  and  partition  function  values, 
therefore  resulting  in  a  more  accurate  vibrational  entropy  contribution.  The 
contribution(s)  from  this  mode(s)  are  then  added  back  in  as  free  internal 
rotations.  Testing  of  this  method  involved  twenty  molecules:  four  without 
internal  rotations  and  16  with  free  or  partially  hindered  internal  rotations. 


CHAPTER  2 


CALCULATIONS 


BASIC  METHOD 

First,  MNDO  was  used  to  calculate  the  molecule's  optimized  geometry  and 
thermodynamic  properties.  The  only  thermodynamic  values  that  were  studied 
were  those  computed  at  298K.  This  was  accomplished  on  all  of  the  test 
molecules.  MNDO  was  run  on  the  VAX  11/780  computers  at  the  United  States  Air 
Force  Academy. 

CORRECTIONS 

The  entropy  values  for  each  of  the  16  molecules  with  free  or  partially 
hindered  internal  rotations  were  thtti  modified  by  dividing  the  MNDO  calculated 
vibrational  partition  value  by  the  spurious  vibrational  mode  contributions  to 
the  partition  functon.  The  spurious  vibrational  mode  contributions  to  the 
vibrational  entropy  were  then  removed  from  the  entropy  value.  Finally,  the 
spurious  contributions  to  the  partition  and  entropy  values  were  replaced  with 
the  corresponding  completely  free  internal  rotation  values. 

The  spurious  vibrational  contributions  to  the  partition  function  were 
computed  using  the  following  equation: 

qvnth  mode  ,  x  /  (1.ev/kT)  (i)3 

where  k  *  6.95196xl0~*  cm”*K  *  and  v  is  the  frequency  of  the  nth  mode 
in  cm  *  (n  ranging  from  1  to  the  number  of  free  IRs).  The  vibrational 
contribution  to  the  partition  function  was  then  corrected  by  dividing  the  MNDO 
vibration  partition  value  by  the  product  of  the  spurious  vibrational 
contributions . 
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The  contribution  to  the  internal  partition  function  by  the  "free"  internal 


rotation  was  then  computed  using  Equation  2: 

Qfree  IR  *  (2.7935/n)(1038IrT> 1/2  <2)4 


with  n  representing  the  symnetry  number  of  the  rotating  group,  I  being  the 
rotational  moment  of  inertia,  and  T  the  temperature  (K).  Following  this 
calculation,  the  maximum  internal  partition  value  was  determined  using  the 
equation 


max 


v*  r 
=  q  x  q  x  Q 


Free  IR 


(3) 


v*  t  ^ 

where  q  is  the  corrected  vibrational  partition  value,  q  is  the 
rotational  partition  value,  and  QFree  IR  is  calculated  from  Equation  2. 

This  value  was  determined  to  be  the  criterion  for  free  internal  rotation 

corrections,  and  its  importance  will  be  discussed  in  the  next  section. 

Next,  the  entropy  value  must  be  corrected  for  the  free  internal 
rotations.  The  general  equation  for  entropy  (equation  4)  results  in  an 

equation  for  the  entropy  for  free  internal  rotations  (equation  5),  which  is 

then  reduced  to  equation  6. 


S  *  n  R  fln(q)  +  T(dln(q)/dT) J 


(4) 
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SFree  IR  "  “  R  + 


T[n/2.7935)(l/1038I  T) 1/2(2.2935/n)(1038I  ) 1/2(l/2)( 1/T1/2) 1]  (5) 


SFree  IR  "  n  R  +  °*5)  (6) 
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In  these  equations,  S  is  the  entropy,  n  is  the  number  of  moles  (assumed  to 
be  one  for  our  purposes),  R  is  the  ideal  gas  constant,  T  is  the  temperature 
(K),  and  q  is  the  partition  function  value  calculated  in  equation  2.  Equation 
7  is  then  used  for  calculating  the  specific  vibrational  mode  entropy 
contributions  and  it  is  derived  from  the  general  equation  for  entropy. 

8„th  node  *  * 

Il(1.43»8.ilth  »  ,-1-‘388''»th/T)/(Tli_e-1.4388voth/T))  (7) 

The  vibrational  contribution  was  then  corrected  by  subtracting  out  the 
spurious  vibrational  contributions  to  the  vibrational  entropy.  The  assumed 
free  internal  rotation  entropy  contribution  was  then  added  back  in  when 
determining  the  corrected  total  entropy  value. 
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CHAPTER  3 


RESULTS 

In  order  Co  test  the  results  of  these  calculations,  the  value  of  */-  1.0 
entropy  units  (1.0  cal/K)  within  the  experimental  value  was  determined  to  be 
acceptable.  Table  I  shows  the  experimental  and  MNDO  entropy  values  for  all 
test  molecules.  It  was  expected  that  the  MNDO  entropy  values  for  molecules 
which  did  not  have  any  free  or  hindered  internal  rotations  would  be  fairly 
close  to  the  experimental  value,  and  all  four  test  molecules  that  met  this 
criteria  (methane,  ethene,  ethyne,  and  propadiene)  were  within  acceptable 
limits.  However,  only  six  of  the  remaining  sixteen  molecules  with  free  or 
hindered  internal  rotations  were  considered  acceptable.  Five  of  these  the  six 
molecules  had  free  internal  rotations  that  were  adjacent  to  a  C-C  single 
bond.  Therefore,  whenever  the  free  internal  rotor  is  adjacent  to  a  double  or 
triple  bond,  MNDO  calculations  are  only  within  acceptable  limits  in  one  out  of 
ten  molecules.  Thus  it  was  concluded  that  hindered  internal  rotations,  such 
as  those  in  saturated  hydrocarbons,  can  be  modeled  as  low  frequency  vibrations 
with  a  small  degree  of  inaccuracy,  while  relatively  free  internal  rotors 
cannot. 


TABLE  I 


MNDO  Entropy  Calculations 


Molecule 

S®xP(cal/K)a 

SM*®°(cal/K) 

Error  (cal/K) 

Methane 

44.50 

44.51 

-0.01 

Ethene 

52.44 

52.29 

+0.15 

Ethyne 

47.99 

47.05 

♦0.94 

Propad iene 

58.30 

57.67 

+0.63 

Ethane 

54.85 

55.20 

-0.35 

Propane 

64.51 

65.45 

-0.94 

Propene 

68.80 

65.45 

+3.35 

Propyne 

59.30 

58.02 

•  +1.28 

n-Butane 

74.10 

73.44 

+0.66 

1-Butyne 

69.51  ' 

74.57 

-5.06 

2-Butyne 

67.71 

73.25 

-5.54 

1-Butene 

73.48 

76.07 

-2.59 

cis-2-Butene 

71.90 

79.52 

-7.62 

trans-2-Butene 

70.86 

87.09 

-16.23 

2 -me thy 1-propane 

70.42 

71.37 

-0.95 

2-methy 1 -propene 

70.17 

80.47 

-10.30 

1,2  Butadiene 

70.03 

70.15 

-0.48 

1,3  Butadiene 

66.62 

68.85 

-2.23 

Methanol 

56.80 

57.05 

-0.25 

Ethanal 

63.50 

66.63 

-3.13 

aAll  experimental  entropy  values  are  from  the  CRC  Handbook  of 
Chemistry  and  Physics,  48th  Ed.,  (Cleveland,  Ohio:  Chemical  Rubber  Co., 
1967),  p.  D-51. 


The  next  important  step  was  to  develop  a  criterion  which  would  determine 
when  the  correction  for  free  internal  rotations  should  be  used.  Since  the 
problem  stemmed  from  internal  rotational  barriers,  an  attempt  was  made  to  find 
the  solution  within  these  barriers  by  comparing  their  values  with  the  errors 
made.  However,  after  compiling  the  necessary  data,  it  was  discovered  that  the 
internal  rotational  barriers  were  very  difficult  to  calculate  or  find 
experimental  data  for,  sometimes  poor  in  agreement  between  experimental 
methods,  and  poor  in  correlation  to  the  results  (see  Table  II)  Although  there 
was  a  general  rule  for  when  to  correct  or  not  to  correct  using  internal 
rotational  barriers  (at  about  2.0  kcal/mol),  this  method  was  discarded  because 
of  its  difficulty  and  lack  of  accurate  results. 


7 


TABLE  II 


Entropy  Differences  and  Internal  Rotation  Barriers 


Molecule 

gCORR 

Error***®® 

Error®0*8 

I.R.  Barrier(kcal/mol) 

Ethane 

55.66 

-0.35 

-0.81 

2.93a 

Propane 

66.40 

-0.94 

-1.90 

3.57a 

Propene 

64.63 

♦  3.35 

+4.17 

1.95a 

n-Butane 

79.88 

+1.28 

-5.78 

3. 25-3. 40 ®( methyl ) 

1-Butene 

77.27 

-2.59 

3.79 

3.16-3.99* 

cis-2-Butene 

76.69 

-7.62 

-4.79 

0.45a 

trans-2-Butene 

72.58 

-16.23 

-1.72 

1.95® 

2 -me  thy 1-propane 

74.75 

-0.95 

-4.33 

3.6-3. 9a 

2-me  thy 1-propene 

71.59 

-10.30 

-1.42  * 

2.12-2.35® 

1,2  Butadiene 

69.94 

•  -0.48 

+0.09 

1.59® 

1,3  Butadiene 

71.25 

-2.23 

-4.63 

2. 3-5.0® 

Methanol 

59.02 

-0.25 

-2.22 

1.07-1. 60b 

Ethanal 

63.54 

-3.13 

-0.04 

1.0-1. 15b 

®Values  found  in  Internal  Rotation  in  Molecules,  Ed.  W.J. 
Orville-Thomas  (New  York:  John  Wiley  and  Sons, 1974),  p.  413-415. 

bValues  found  in  Conformational  Analysis,  Ernest  L.  Eliel,  (Washington 
D.C.:  American  Chemical  Society,  1981),  p.  140. 


The  next  method  studied  concerned  the  ratio  between  q1  caicuiate<j  and 

q1  .  (See  Table  III).  If  the  program  describes  the  internal 

^  maximum 

rotation  incorrectly  as  a  vibrational  frequency  which  is  too  low,  then  it  will 
show  up  as  a  ratio  that  is  grater  than  one.  If  the  internal  rotation  can  be 
modeled  as  a  low  vibrational  frequency,  such  as  a  hindered  internal  rotation 
can,  then  it  will  show  up  as  a  ratio  that  is  less  than  one,  indicating  that 
the  MNDO  calculation  is  within  the  maximum  bonds. 

This  method  worked  correctly  for  13  of  the  16  test  molecules,  incorrectly 
for  only  two,  and  didn't  matter  for  one  molecule  (see  Table  IV).  In  addition 
to  the  high  percentage  of  correct  decisions,  the  average  percent  deviation 
from  experimental  entropy  values  using  this  criteria  is  only  2.52Z,  while 
using  the  correction  at  all  times  is  3.92Z  and  the  uncorrected  MNDO  results 
give  a  5.50Z  deviation. 
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TABLE  III 


Internal  Partition  Function  Ratios 


Molecule 

9l  MNDO 

91  max 

9*  ratio 

Correction(A/T) 

Ethane 

1.4486x103 

3.1651x10s 

0.4577 

0.46  A 

Propane 

4.5528xl04 

1.7666xl0S 

0.2577 

0.96  A 

Propene 

7.3482x1 04 

6.9547X104 

1.0566 

0.82  A 

Propyne 

3.7918x10" 

1.2126x10; 

0.3127 

0.10  A 

n-Butane 

4.4937x10s 

2.1572x10' 

0.0208 

4.50  A 

1-Butyne 

2.2538xl06 

7.1291x10® 

0.3161 

0.00  A/T 

2-Butyne 

7.6486x10s 

1.2377x10s 

6.1797 

4.66  T 

1 -Butene 

2.8119x10® 

9.5216xlOP 

0.2953 

1.20  A 

cis-2-Butene 

1.1264xl07 

6. 1673x10® 

1.8264 

2.83  T 

trans-2-Butene 

3.6809x10® 

6.4725xlOS 

568.6980 

14.51  T 

2 -me  thy 1-propane 

1.7454x10s 

1.5044x10® 

0.1160 

3.38  A 

2 -me thy 1-pr opene 

1.9828xl07 

5.1678x10s 

38.3684 

8.88  T 

1,2  Butadiene 

2.7893x10s 

2.9513x10s 

0.9451 

0.39  T 

1,3  Butadiene 

1.5139x10s 

7.6343x10s 

0.1983 

2.40  T 

Methanol 

4.3509x10s 

1.1642X104 

0.3737 

1.97  A 

Ethanal 

1.5498x10s 

5. 1427x10^ 

3.0136 

3.09  T 

CHAPTER  4 


RECOMMENDATIONS 


There  is  still  more  to  be  accomplished  before  these  results  can  be 
incorporated  into  a  program  such  as  MNDO.  Since  most  of  the  testing  was  done 
with  hydrocarbons  and  with  free  internal  rotations  which  were  predominantly 
methyl  groups,  it  is  important  to  verify  this  test  with  different  internal 
rotors,  such  as  NH2,  halogen,  and  deuterium  substituted  methyl  groups,  etc. 
If  more  tests  conclude  that  this  correction  is  accurate  and  reliable  enough, 
then  the  results  can  be  incorporated  into  MNDO.  Currently,  a  stand-alone 
version  of  the  subroutine  THERMO  from  MNDO  with  these  corrections  in  it  is 
being  used  to  determine  whether  the  separate  internal  rotation  corrections 
for  the  vibrational  contribution  to  entropy  should  be  used  or  not. 
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